Stability against a decay of some recently observed superheavy elements 
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The probability of a particle emission for some recently observed superheavy nuclei (SHN) are 
investigated. The a-decay half lives of SHN are calculated in a quantum tunneling model with 
density dependent M3Y (DDM3Y) effective nuclear interaction using theoretical and measured Qa 
values. We determine the density distribution of a and daughter nuclei from the relativistic mean 
field theory (RMF) using FSUGold force, NL3 and TMl parameter sets. The double folded nuclear 
potential is numerically calculated in a more microscopic manner using these density distributions. 
The estimated values of Q-decay half-lives are in good agreement with the recent data. We compare 
our results with recently detected a-decay chains from new element with atomic number Z=117 
reported by JINR, Dubna. Finally, we determine the half-lives of superheavy elements with Z=108- 
120 and neutron number N= 152- 190 to explore the long-standing predictions on the existence of an 
"island of stability" due to possible spherical proton (Z~ 114) and neutron (N~ 184) shell closures. 

PACS numbers: 27.90. -(-b, 23.60.-|-e, 24.10.Jv, 21.30.Fe 



The formation of superheavy nuclei (SHN) by fusion 
was intensively explored [IHll- The investigation of ex- 
perimental data concerning fusion and fission of the SHN 
with Z — 108-118 together with data on survival prob- 
ability of these nuclei in evaporation channels with 3-4 
neutrons, revealed the fact that the hindrance due to 
high fission barriers causes to a relatively higher stability 
of such heavy nuclear systems [3]. The investigation of 
a-decay chains of such SHN is the main tool to extract 
some information regarding their degree of stability and 
possible existence in nature. In the last decade several 
theoretical and experimental works were devoted to 
the formation of SHN and their a-decay half-lives. The 
a decay of superheavy nuclei possible if the shell 

effect supplies the extra binding energy and increases the 
barrier height of fission. 

In our previous works [1, [l^, [l^ we showed the ap- 
plicability of our calculation using DDM3Y interaction 
in predicting the a decay half lives of SHN from a di- 
rect comparison with the experimental data 0, 
this paper, we calculate the a-decay half-lives of new su- 
perheavy element Z=117 and its decay products recently 
observed in JINR, Dubna [l7|- The density distribution 
of a-particle and daughter nuclei are determined from 
the RMF theory using FSUGold force [3, NL3 and 
TMl [20| parameter sets. The Lagrangian model NL- 
SVl with the inclusion of the vector self-coupling of the 
w-meson was successfully used by Saldanha, Farhan and 
Sharma [2l'| to study the properties and shell structure 
of the superheavy elements from Z=102-120 based on 
RMF-I-BCS calculations for an axially deformed config- 
urations of nuclei. To include the deformation effect to 
some extent, we have used the Qa values obtained from 
the force NL-SVl [l^l for the investigation of existence 
of the so-called "magic island" with extra stability. How- 
ever, we have determined density distribution from RMF 



calculation assuming spherical nuclei. In addition, the 
expression of centrifugal barrier is modified to introduce 
an additional turning point near the origin, even for the 
I = case. 

The half lives for a disintegration process [23j are cal- 
culated using WKB approximation [2^ in the frame work 
of quantum mechanical tunneling of an a particle from 
a parent nucleus [l^ . The details of calculation of the a 
decay half lives (Ta) of superheavy nuclei were described 
in our earlier works [l^ . The required nuclear poten- 
tials are calculated by double folding the nuclear density 
distribution functions of the a particle and the daughter 
nucleus with DDM3Y effective interaction. The values of 
parameters used in the DDM3Y interaction are k ept un- 
changed as calculated in our previous papers (TBI . [25| . In 
this work, the density distributions are obtained from the 
nuclear wave functions (tp) using the Lagrangian param- 
eter sets FSUGold, NL3 and TML The microscopic a- 
nucleus potential thus obtained, along with the Coulomb 
interaction potential and the minimum centrifugal bar- 
rier required for the spin-parity conservation, form the 
potential barrier. The spin-parity conservation forces a 
minimum angular momentum (/) to be carried away in 
the decay process. Consequently, this gives rise to the 
Langer modified centrifugal barrier 
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where /i is the reduced mass of the daughter and emitted 
nuclei system and R is the distance between them. In 
this work, minimum centrifugal barrier is assumed by 
using I — 0. The Langer modification from l{l + 1) ^ 
{I + 1/2)^ is a necessary transformation while going 
from the one dimensional problem with x ranging from 
— oo — !■ +00 to the radial one-dimensional tunneling with 
r ranging from — >■ c». 

In this calculations, the experimental Qa are used for 
the estimate of of the recently synthesized new ele- 
ment Z=117 in JINR, Dubna [l3l and recently observed 
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decay of 288,289 ;lj^4 Qgj^ Darmstadt [13 • The the- 
oretical Qc-values for the neutron rich isotopes of the 
elements with even atomic number Z=108-120 are ob- 
tained from the force NL-SVl to search the existence of 
an "island of stability" surrounded by unstable nuclear 
systems. 

The experimental decay Q values {Qex) have been ob- 
tained from the measured a particle kinetic energies Ea 
using the following expression 



)Ea + (65.3Z^/^ - 80.0^2/^) X 10"^ MeV 

(2) 

where the first term is the standard recoil correction and 
the second term is an electron shielding correction. The 
mass and atomic numbers of the parent nucleus are de- 
noted as Ap and Zp respectively. As the -value ap- 
pears inside the exponential integral as well as in denom- 
inator of the expression [l^ of a decay half lives, the 
entire calculation is very sensitive to Q- values. 

The RMF theory has successfully described various 
properties of nuclei frorn light to superheavy domain. 
The RMF calculations are used to explore the na- 
ture of possible magic numbers using various interactions 
and it is shown in some calculations that the patterns 
of single -particle levels are significantly modified in su- 
perheavy elements [l^. Being based on the Dirac La- 
grangian density (£), RMF is particularly suited to in- 
vestigating these nuclei because it naturally incorporates 
the spin degrees of freedom. There are different varia- 
tions of the Lagrangian density and also a number of dif- 
ferent parameterizations in RMF namely, FSUGold force, 
NL3 and TMl etc. The following Lagrangian density has 
recently been proposed |18j which involves self-coupling 
of the vector-isoscalar meson as well as coupling between 
the vector-isoscalar meson and the vector-isovector me- 
son. The corresponding parameter set is called FSUGold. 

C = ^P (Z7^a'^ - Mn) V- + ^ {d^ad^'a - mla^) 
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d^A, - d^A^ (4) 
- a^w,, - duUJ^ (5) 
dvPt, - Qp {PpXp^) (6) 
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In this work, we determine the density distributions of 
emitter and daughter nuclei from RMF calculations us- 
ing three different parameter sets namely, FSUGold force 
m. NL3 [H and TMl [H. 



In this short report, our primary aim is to check the 
stability of recently discovered jlTI] new superheavy ele- 
ment with atomic number Z=117 against a-decay. The 
two new isotopes ^'^^117 and ^^'^117 were recently pro- 
duced along with 9 more new nuclei as their decay prod- 
ucts at JINR, Dubna in fusion reaction between ''^Ca 
and 249gjj_ |.]-^jg work, the a-decay lifetimes (r) are de- 
termined within a WKB framework using density distri- 
butions for daughter and emitted nuclei from RMF calcu- 
lation. In Table I the experimental {t^-^^) [i3] and the- 
oretical a-decay lifetimes {Tp^suGoicn ^nl3 ^ '^tmi') calcu- 
lated by using density distributions from FSUGold, NL3 
and TMl parameter sets are compared for 293,294-[^-|^y ^^^^ 
their decay products. In addition, the calculations are 
also done with theoretical Qa values used in Ref. [13] 
based on macroscopic-microscopic calculations of masses 
of superheavy nuclei [soj . The upper and lower limits 
of calculated values are arisen due to the error encoun- 
tered in measured Qa values. It is clear from Table I that 
the three sets of calculated a-decay hfctimcs Tp^suGoid^ 
TnlI^ Ttmi using wave functions from FSUGold, NL3 
and TMl parameter sets respectively with DDM3Y in- 
teraction are close to each other. Our calculations using 
experimental Qa are in reasonable agreement with exper- 
iments for some nuclei like "^^^117, ^s^llS, '^^'^117, 2*^113 
etc. For the rest of the nuclei like ^soilS, 282iii^ 278^09^ 
2 ''^107 etc. our calculations are in better agreement with 
theoretical lifetimes given in Ref. [iTj . For example, the 
calculated lifetimes using FSUGold 11.3 seconds) is 
much less than measured lifetimes 1.3 minutes) for 
2^^107. However, this discrepancy is removed if we em- 
ploy the Qa values from Ref.[30|. The data with higher 
statistics are needed to determine the lifetimes of SHN 
with a better accuracy. 

Recently the a decay chains from 289214 and 288j^j^4 
are reported by DuUmann et at. [13]. In Table II the 
calculated a decay half lives {T[^l^°°^'^, T^^^^, T™^) 

and experimental a decay half lives (Tjy^^) are given 
for comparison. The half lives predicted in this work are 
in good agreement with the measured values for 289 114 
and 2*i£)s. For nuclei 285ii2 and "^^^114 our calculations 
are in reasonable agreement with experiments. 

Finally, the a-decay stability of a series of superheavy 
nuclei with Z=108-120 towards neutron rich domain are 
investigated by using Qa values obtained from the force 
NL-SVl [2l[. The nuclear potential are numerically de- 
termined using double folding the DDM3Y interaction 
with density distribution functions obtained from the 
FSUGold parameter set. The calculated halfiives are 
shown in Fig. 1 with increasing neutron number N=152- 
190 for elements with even atomic number Z=108-114. 
The half life for any specific element (i.e. for fixed Z) 
increases with increasing neutron number. This matches 
the trend of higher stability towards neutron shell closure 
around N~ 184 or 186 as predicted by many microscopic 
calculations. A sharp fall in Ta around N=184 for ele- 
ments Z=108, 110, 112, 114 may indicate the signature of 
shell effects around N=184 or 186. A similar trend is ob- 
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set. It seems from this calculation that the longer lived 
superheavy nuclei may exist due to shell effect near 
184 with Z=ll 0-114. However, more experimental 
and theoretical works are needed to be done to shed 
more light on this issue. 
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FIG. 1. (Color online) Half lives log T^, (in second) vs. neu- 
tron number (N) using density distribution from the force 
FSUGold and Q„ from NL-SVl parameter set for Z=108-114 



The research work of P. Roy Chowdhury is sponsored 
by the UGC (No.F.4-2/2006(BSR)/13-224/2008(BSR)) 
under Dr. D.S. Kothari Postdoctoral Fellowship Scheme. 
The author (P.R.C.) is thankful to Prof. M.M. Sharma 



4 

2 

-2 
-4 - 
-6 - 
-8 
-10 




150 



160 



170 180 
N 



190 



200 



served in Fig. 2 for elements Z=116, 118, 120. However, 
as the atomic number increases the shell effect seems to 
be less prominent. For example, the shell effect near 
N=184 becomes less significant for Z=118 and 120. 

In summary, the a-decay halflives Tq, of some recently 
observed superheavy nuclei are determined in a more 
microscopic approach using DDM3Y interaction and 
RMF wave functions obtained from three different pa- 
rameter sets FSUGold, NL3 and TMl. The calculations 
are in reasonable agreement with the measured values 
of Tq,. Finally, a thorough investigations are done to 
search the existence of an "island of stability" in neutron 
rich domain of superheavy nuclei with Z=108-120 and 
N=152-190 using the Qa values from NL-SVl parameter 



FIG. 2. (Color online) Half lives log (in second) vs. neu- 
tron number (N) using density distribution from the force 
FSUGold and from NL-SVl parameter set for Z=l 16-120 
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TABLE I. Comparisons between observed (r^"^ ^) [V^ and theoretical (this work) a-decay lifetimes (t^s^Jq^j^, t^'U , '''tmi) 



using measured and calculated Qa given in Ref.[17l]. The calculated Qa values taken from Ref. [3C| are given with single (*) 
asterisk symbol. The respective half lives given in Ref.[l7| using calculated Qa are also denoted by single (*) asterisk symbol. 

Parent Ref. fl7] Ref. fl7] Lifetimes Lifetimes Lifetimes Lifetimes 

E-''* (MeV) g. {MeV) r^^^ Ref.fl7] T^MZ, ^ 



293^17 


11 03 f8~l 


11.19 (8) 


21 ms 


4.1 + 2-3 jy^g 


3.5+2-1 


3.6+2-1 


293^17 


*11.26 


11.42 


*10 ms 


1.2 ms 


1.0 ms 


1.0 ms 


28^115 


10 31 fq"! 


10 46 fQ") 


0.32 s 


65+*7 ™„ 


56+*" ms 


57+41 


289115 


* 10.48 


10.63 


*0.22 s 


23 ms 


20 ms 


20 ms 


285 


9.48 (11) 


9 62 n 1 "1 


7.9 s 


3 0+^* s 

'^•'-'-1.6 * 


2 6+2-8 s 

— 1.4 


2 6+^2 „ 

— 1.4 


285113 


9.74 (8) 


9.89 (8) 


7.9 s 


530t^^^ ms 


456j:?82 


464+?" ms 

— 187 


285 


*9.96 


10.11 


*1.2 s 


0.13 s 


0.11 s 


0.11 s 


29*117 


10.81 (10) 


10.97 (10) 


112 ms 


521*1 ms 


45l2o 


48^22 


ii ( 


11. UU 


11. ID 


*45 ms 


lo TflS 


15 ms 


16 ms 


290115 


9.95 (40) 


10.10 (10) 


0.023 s 




2.09lf,« s 


9 -,0+28.89 
^•^•^-1-96 * 


290115 


*10.23 


10.38 


*1.0 s 


0.40 s 


0.34 s 


0.35 s 


286113 


9.63 (10) 


9.77 (10) 


28.3 s 


4.2l*;0 s 


3.61?:* s 


3.6i?;^ s 


286113 


*9.56 


9.70 


*16 s 


6.7 s 


5.7 s 


5.8 s 


282111 


9.00 (10) 


9.13 (10) 


0.74 s 


70 1+77 3 s 


59.95I3Q yg S 


60.9lt«I:«? ^ 


282111 


*9.43 


9.57 


*8.1 s 


3.4 S 


2.9 S 


3.0 s 


2™ 109 


9.55 (19) 


9.70 (19) 


11.0 s 


0.3tO:« s 


o.3«:« s 


o.3t«:^ s 


278 -^Qg 


*9.14 


9.28 


*13 s 


5.0 s 


4.3 s 


4.3 s 


27*107 


8.80 (10) 


8.94 (10) 


1.3 min 


11.3tjy s 




9.7t^°o' s 


27*107 


*8.43 


8.56 


*7.4 min 


3.02 min 


2.57 min 


2.62 mm 



TABLE 11. Comparisons between observed ^ and theoretical (this work) a-decay half lives (^rf^fuGoid^ , T™^) using 

measured Qa [27| . 



Parent 


Ref. [27] 


Ref. [27] 


Halflives 


Halflives 


Halflives 


Halflives 




E:^^p' (MeV) 


Qa [MeV) 


T^/V Ref. [27] 


rpFSUGold 

-'1/2 


rpJVL3 

-'1/2 


rrnTAfl 
-'1/2 


289114 
285112 

28lDs 
288114 


9.87 (3) 
9.21 (3) 
8.73 (3) 
9.95 (3) 


10.06 (3) 
9.39 (3) 
8.90 (3) 
10.14 (3) 


0.97l°;^7 s 
30l?o s 
140l^jo s 
0.471°;?* s 


n or+0-08 

6.40l?:fo ^ 
Ai.&ltf.^J s 

0.091^2 - 


n OO+0-07 

6.05ll:« s 
39.18+?-^? s 

o.ostH? - 


31+0 07 

^•■31-0-05 ^ 

5.65l?:gi s 
38.08+7 g s 
0.0810:01 s 
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